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We study a spatially extended model of noise-induced 
magnetization reversal: a classical Ginzburg-Landau the- 
ory, restricted to a bounded interval and perturbed by weak 
spatiotemporal noise. We compute the activation barrier and 
Kramers prefactor. As the interval length increases, a tran- 
sition between activation regimes occurs, at which the pre- 
factor diverges. We relate this to transitions that occur in 
low-temperature quantum field theory. 

PACS numbers: 05.40.-a, 05.45.Yv, ll.10.Wx, 75.60.Jk 



The effect of noise on spatially extended classical- 
mechanical systems has become a subject of intense in- 
vestigation Q. Noise, by which is meant local fluctua- 
tions of thermal or other origin, may be spatiotemporal: 
it may vary randomly in space as well as time. 

A typical problem is the determination of the extent 
to which spatiotemporal noise can induce transitions be- 
tween the stable states of a system described by a non- 
linear field equation, or induce such a system to escape 
from a metastable state. Langer || and Biittiker and 
Landauer j^] studied one-dimensional versions of this 
problem, in the limit of an infinite-size spatial region. 
Recently, finite size effects have become a subject of 
study and simulation. Micromagnetics is one applica- 
tion area. The thermally activated magnetization rever- 
sal of nanoscale magnets is described by coupled Landau- 
Lifschitz-Gilbert equations perturbed by spatiotemporal 
noise |Q. An interesting question is how the effects of 
spatial extent cause magnetization reversal to differ from 
the better-understood 'zero-dimensional' case of a single- 
domain particle, with the noise taken to have no spatial 
dependence j^] . Another area where activation by spatio- 
temporal noise is important is the formation of spatially 
localized structures in electroconvection ||. 

In this Letter, we quantify the effects of weak spatio- 
temporal noise on an overdamped bistable quartic (i.e., 
double well) classical Ginzburg-Landau field theory, in 
a one-dimensional region: the bounded interval [0, L]. 
This model is similar to that of Langer Q , and has much 
in common with the sine-Gordon model of Biittiker and 
Landauer S. For sufficiently large L, the model has 
two stable states, which are states of positive and neg- 
ative magnetization. In the Kramers weak- noise limit, 
in which noise-activated magnetization reversals become 
exponentially rare, we compute the reversal rate T. We 
also determine the 'optimal trajectories', in the model's 
infinite-dimensional state space, for magnetization rever- 
sal. They are greatly affected by the choice of boundary 



conditions. If periodic boundary conditions are used, it is 
most likely that reversal will proceed via the nucleation 
of a droplet within a pair of Bloch walls; but for Dirichlet 
or Neumann boundary conditions, a single wall will form 
at x — or x — L, and sweep across the interval. Faris 
and Jona-Lasinio worked out a mathematically rigorous 
'large deviation theory' of the Dirichlet-case reversal M. 

Reversals fall off according to T ~ To exp(— AE'/e), 
where e is the noise strength, AE is the activation barrier 
between the two states, and To is the Kramers prefactor. 
We show how to compute the prefactor in closed form, 
as a function of L, for all three boundary conditions, 
and in the Dirichlet case, work it out in full. Langer 
computed it in the L — > oo limit of the periodic case, 
so this is a significant advance. The prefactor contains a 
quotient of infinite-dimensional fluctuation determinants. 
Since the work of Coleman pj , the quantum field theory 
community has known how to compute such quotients as 
L — > oo. With the regularization technique of McKane 
and Tarlie M, we can handle the L < oo case, as well. 

A feature of our prefactor computation is the calcu- 
lation of the unstable eigenvalue of the model's deter- 
ministic dynamics, at the transition state between its 
two stable states. The Kramers prefactor differs sig- 
nificantly from the analogous low-temperature quantum- 
mechanical tunneling prefactor, in that it involves this 
eigenvalue jlO| . Calculating it in the Dirichlet case is dif- 
ficult. We express it in terms of a mid-band eigenvalue of 
the / = 2 Lame Hamiltonian fll]Jl2]| . Using a dispersion 
relation developed for this Hamiltonian |iljf| , we compute 
the eigenvalue, and hence the Dirichlet-case prefactor. 

Our most striking discovery has to do with the phase 
structure of the weak-noise limit. In a sense that can be 
made precise, the e — > limit of the quartic Ginzburg- 
Landau theory has a second-order phase transition at 
L = 2tt (Dirichlet, periodic cases) or L = it (Neumann 
case), in dimcnsionless units. At these critical values 
of the interval length, the Kramers prefactor diverges. 
This is due to the bifurcation of the transition state, as 
L increases through criticality. A zero-field 'sphaleron' 
configuration, which serves as transition state when L is 
small, bifurcates into a degenerate pair of 'periodic in- 
stantons'. Similar bifurcations have been studied by the 
quantum field theory community |l4]-[l7|] . Until now, the 
importance of these bifurcations to classical activation by 
spatiotemporal noise has not been fully realized. 

Model and Phenomenology. — In specifying the over- 
damped Ginzburg-Landau model, we follow Faris and 
Jona-Lasinio Q. On [0, L], a classical field <j> — <p{x,t) is 
evolved by the stochastic Ginzburg-Landau equation 
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M cj>" + ^ - \(f) 3 + e 1/2 £(x, i), 



(1) 



where £(x,t) is unit-strength spatiotemporal white noise, 
satisfying {£(xi, h)£(x 2 , t 2 )) = S(%i - x 2 )S(t 1 - t 2 ). 
We set M = fj, = A = 1, or equivalently use dimensionless 
units, in which the length unit is the nominal coherence 
length \J Mj [i, the field strength unit is y^^/X, etc. 

In the absence of noise, the time-independent solutions 
of (Q) include the sphaleron <p = 0, and cf> = ±1; the latter 
are the magnetized states in the Neumann and periodic 
cases. In the Dirichlet case, 4>{x = 0) = 4>(x = L) = 
is imposed, so the magnetized states necessarily differ. 
It is easy to check that if e = 0, the so-called periodic 
instanton <p = 4>inst,m(x) is a time-independent solution 
of (0) for any to in the range < m < 1. Here (cf. 
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i(x/ vm+T | to), 



(2) 



where sn(» | to) is the Jacobi elliptic function with pa- 
rameter to, which has quarter-period K(m) [jl8j. This 
quarter-period decreases to ir/2 as to — > + , in which 
limit sn(» | to) degenerates to sin(»), and it increases to 
infinity asm^l". One would expect the Dirichlet-case 
stable magnetized states to be ±<?!>inst,m,, D , with m s> D de- 
termined implicitly by the condition <p(x = L) = 0, i.e., 
by the half-period condition 2y/m Si rj + lK(m Sj fl) = L. 
However, this is so only if L > ir. If L < ir, there is 
no solution for to Sj £> in the range < m s u < 1, and 
the Dirichlet-case model is monostable with the zero- field 
configuration as its stable state, rather than bistable. 
Bistability disappears when L — > ir + and m Si rj — » + . 

Noise- activated magnetization reversal, for weak noise, 
proceeds with high likelihood along an optimal trajectory 
in the model's infinite-dimensional state space that goes 
'uphill' from a magnetized state to an intermediate tran- 
sition state. In a zero-dimensional approximation, this is 
the sphaleron. However, on physical grounds one expects 
a transition state to be an unstable field configuration 
with two droplets, i.e., with half the interval occupied by 
each magnetization value. Mathematically, this can arise 
as follows. In the Dirichlet case, the sphaleron undergoes 
a pitchfork bifurcation into the degenerate magnetized 
states ±(/>i nst . ms D when L is increased through n. When 
L is increased through 2-7T, it bifurcates again, into a de- 
generate pair of unstable states ±<fii nstiTnu D . Unlike the 
magnetized states, these transition states are positive on 
half the interval and negative on the other. m Ut D is de- 
termined from L by the condition <p(x = L) = 0, i.e., by 
the quarter-period condition 4^/r7i u ^ + 1 K(to u .£>) = L. 

The same bifurcation occurs if periodic boundary con- 
ditions are used, but the new transition state is infinitely 
degenerate: it is 0inst,m„ D , translated arbitrarily. It is 
easy to see that in the Neumann case, the bifurcation 
occurs at L = n rather than L = 2tt, with the Neumann 
parameter m u ^ determined by the half-period condi- 
tion 2y/m u ^ + 1 K(m n .Ar) = L. In fact, m u ,N — m-s.D- 
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FIG. 1. Stable magnetized states (S) and transition 
states (U), for Dirichlet (D), Neumann (N), and periodic (P) 
boundary conditions, if L = 10. Each state may be multiplied 
by —1, and UP may be shifted arbitrarily. 

To satisfy Neumann boundary conditions, the field con- 
figurations ±<j)- lnst .m u jv must be shifted by L/2. 

Figure 1 displays the magnetization and transition 
states for each boundary condition. As L — > oo, the 
Bloch walls between magnetization values increasingly 
acquire the standard hyperbolic tangent form, since 
sn(» | to) degenerates to tanh(«) as to — > 1~. 

The zero-noise dynamics of this model are of the gra- 
dient form 4> — —SH/S<p, with the energy functional 



[(0') 2 /2 + ^ 4 /4 - 2 /2 + 1/4] dx. (3) 



The energy of each magnetization and transition state is 
readily computed from (||). The states <f> = ±1 have zero 
energy, and the sphaleron has energy L/A. The energy 
of each periodic instanton state ±</>j ns t,m turns out to be 
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8 E(m) (1-to)(3to + 5) 



(m + 1) K(to) 



(to + l) 2 



(4) 



where E(m) is the second complete elliptic integral |lq| . 

Figure 2 plots the activation barrier AE as a function 
of L, for each boundary condition. The second deriva- 
tive of each AE function is discontinuous at the value 
of L at which bifurcation occurs. As L — > oo, the value 
of AE converges to the energy of a Bloch wall (Dirich- 
let, Neumann cases), or two Bloch walls (periodic case). 
It follows from (||) that the energy of a Bloch wall, of the 
to -> 1" limiting form tanh(x/>/l), is 2^2/3 w 0.943. 

Determinant Quotients. — The formula for the Kramers 
prefactor Tq of an overdamped multidimensional system 
driven by weak white noise is well known. Suppose the 
system has a stable state ip s and a transition state <p u , 
with a single unstable direction. Let A s and A u denote 
the system's linearized noiseless dynamics at <p s and ip u , 
so that to leading order, the state ip = ip s + 1] evolves by 
f\ = A s r], and ip = tp u + r\ by 77 = -A u r]. Then Jl9| 



2 



2.0 



1.5 



AE 1.0 



0.5 



0.0 









N 











10 



FIG. 2. The activation barrier AJ5, for the cases of Dirich- 
let (D), Neumann (N), and periodic (P) boundary conditions. 
Bullets indicate criticality (L = 2n for D and P; L = -n for N). 
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T| |A„,i| 



det A, 



det A u 



(5) 



where X Ut % is the only negative eigenvalue of A„. The cor- 
responding eigenvector rj u 1 is the direction along which 
the optimal trajectory approaches the transition state. 
If A Sj i denotes the smallest eigenvalue of A s , the cor- 
responding eigenvector 77 s 1 will be the direction along 
which the trajectory extends from the stable state. 

Linearizing the noiseless version of (]]]) at a stationary 
state <fio (either a stable or a transition state) yields 
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-A[0o] T) 



[-d 2 /dx 2 + (-l + 3(j> 2 Q )]r h (6) 



So Tq depends on the spectrum of the A operators as- 
sociated with the stable and transition states. In the 
Dirichlet case, the formal determinant quotient T can 
be computed by the now standard technique of Cole- 
man (|||. If > 27r, let r? Si * and rj u ^ be the solutions 
on [0, L] of the homogeneous equations A.[4>[ nst nis D ]rj = 
and A[(f>i nstj m u D ]v = which satisfy the boundary con- 
ditions 77(0) = and 7/(0) = 1- Then, it turns out, 
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Va,*{L)/ri Ut *(L). 



(7) 



is the Dirichlet-case determinant quotient. 

Solutions r) s .„ and r) u> * satisfying these special bound- 
ary conditions may be constructed by a clever trick ||: 
differentiating the periodic instanton </>i n st,m with respect 
to to, setting m to m St £> and m Ut £> respectively, and nor- 
malizing. This procedure uses the formula for the deriva- 
tive of sn(» I to) with respect to m |2(|. The result is 



L 



m c n + 1 



.,£) + ! E(m c<D ) 



1 - m c , D K(m Ci rj) 



for both c = s and c — u ('±' is positive and negative 
respectively). Substituting rj s ^(L) and rj u ^(L) into (]?]) 
yields the determinant quotient To- 



The unbifurcated regime, i.e., n < L < 2ir, must be 
handled a bit differently. Since the transition state is 
the sphaleron, and not the periodic instanton <fri ns t,m u ,n > 
Vu,*{L) cannot be computed from the above formula. But 
it is trivial to check that rj u ^(L) simply equals sinL. 

The Neumann case is treated similarly to the Dirichlet 
(we omit the details), but the case of periodic bound- 
ary conditions is very different, at least in the bifurcated 
regime L > 2tx. The periodic instanton transition state 
is infinitely rather than doubly degenerate, and at any 
transition state, the linearized dynamical operator A has 
a soft 'collective mode', with a zero eigenvalue. Equa- 
tion (]^) must be modified, but the regularization tech- 
nique of McKane and Tarlie (|] can be used to work out 
the periodic-case Kramers prefactor [^fj. It acquires an 
e" 1 / 2 factor. As a result, the periodic-case reversal rate 
becomes non-Arrhcnius when L is increased through 2n. 
A similar non-Arrhenius reversal rate falloff is displayed, 
in the limit of zero spatial extent, by the stochastic 
Landau-Lifschitz-Gilbert equation ||, due to its space 
of magnetization values having a continuous symmetry. 

The Unstable Eigenvalue. — The stumbling block in the 
analytic computation of the Dirichlet-case Kramers pre- 
factor is the calculation of A Ul i, the single negative (un- 
stable) eigenvalue of the deterministic dynamics, lin- 
earized at the transition state. If n < L < 2n and the 
transition state is the sphaleron, X Ut \ = ir 2 /L 2 — 1 is triv- 
ial to verify. But computing X Ut i and the corresponding 
eigenfunction t] Ui i when L > 2ir is much harder. t] Ui i is of 
considerable physical interest, since it describes the way 
in which the optimal trajectories approach the periodic 
instanton solutions ±0inst,m„ D , i.e., the way in which the 
moving Bloch wall slows to a halt at x = L/2. 

Here we sketch the calculation of A Ui i and r\ u ^\ from 
the eigenvalue equation A77 = A77; details will appear 
elsewhere Introducing z = xj ^Jm u ^ + 1, and for 

simplicity, writing m for m n .£>, converts the equation to 



[— d 2 /dz 2 + 67nsn 2 (z | m)] 77 = £ ?/, 



(8) 



where the 'energy' £ equals (m + l)(A U) i + 1). The inter- 
val < z < 4K(to) corresponds to [0, L\. Equation (||) is 
the I — 2 Lame equation |l^JT^] ; which is a Schrodinger 
equation with a periodic potential, whose lattice constant 
is 2K(m). Its Bloch wave spectrum is known to con- 
sist of three energy bands p2| , each extending over the 
wavenumber range — 7r/2K(TO) < k < 7t/2K(to). 

According to Hermite's solution of the Lame equa- 
tion B2] , Eq. (H) has solutions of the form 



H(z ± on I m) 



Q(z I m) 



exp[=F Z(a 2 | m) z] \ , (9) 



where H, 0, and Z are the Jacobi eta, theta, and zeta 
functions, and ai,cv2 are certain complex constants de- 
termined by £. Clearly, the wavenumber k of the solu- 
tion (j|) equals ±Im^ 2 =1 Z(cti \ m). 
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into account, in simple models of magnetization rever- 
sal induced by weak noise, may yield a rich structure of 
activation regimes separated by phase transitions. The 
choice of boundary conditions plays a major role. We ex- 
pect more sophisticated models, and the phenomenon of 
field-induced reversal, will display similar structure. 

This research was partially supported by NSF grants 
PHY-9800979 and PHY-0099484. 
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FIG. 3. The Kramers prefactor To in the Dirichlet case. 

The difficulty lies in finding closed-form expressions for 
a\, 0L2 which will yield a solution r\ = T) u i{ z ) that satisfies 
Dirichlet boundary conditions on < z < 4K(to). But 
Tl u i, being a ground state, should have no nodes, so it 
should extend to a periodic function with period 8K(m). 
That is, its wavenumber should be ±7r/4K(m). One of 
us [ p"3[ has found expressions for a%,a2 in terms of £, 
which yield a dispersion relation £ i— > ±fc. By in- 
verting this, we obtain an energy £ corresponding to 
k = ±7r/4K(m), and hence the eigenvalue X u ,i- 

The Prefactor. — The Dirichlct-case Kramers prefac- 
tor To can be computed from the formula (||), using the 
closed- form expression (^) for the determinant quotient, 
and the just-explained technique of calculating the un- 
stable eigenvalue A U) i< Figure 3 shows the dependence 
of To on the interval length L. The divergence at L = 2tt 
is typical of a second-order phase transition. 

The divergence can be viewed as arising from the 
bifurcation of the optimal weak-noise reversal trajec- 
tory, rather than the bifurcation of the transition state 
on which it terminates. We previously found a simi- 
lar prefactor divergence in a symmetric two-dimensional 
nonequilibrium model, whose optimal trajectory bifur- 
cates but whose transition state does not (§3] . 

A natural question is whether the second-order phase 
transition is robust. In slightly more complicated classi- 
cal field theories perturbed by spatiotemporal noise, such 
transitions may be first-order rather than second-order, 
with a discontinuous, rather than diverging, prefactor. 
Kuznetsov and Tinyakov Jl5| have studied stationary 
field configurations in a sixth-degree Ginzburg-Landau 
theory, with (fi+arf) 3 — (a + l)</> 5 replacing the tfi—cj) 3 terms 
of (0). The a = — 1 case is the theory we have treated. 
If a > —1, the periodic instanton branch of the energy 
function crosses the sphaleron branch at a nonzero angle. 
This would give rise to a first-order transition. So, in the 
(L,a) plane, the second-order transition point (2ir, — 1) 
must be the endpoint of a first-order transition curve. 

To sum up, we have shown that taking spatial extent 
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